clc,clear;
x1=1900:1:2300;
x2=1800:-1:1400;
y1=10:0.01:14;
y2=9:-0.01:5;
for i=1:401
    for j=1:401
Y{1}(i,j)=Theorem3(1817.2,0.53512,0.30758,1800,x1(i),1800-x2(j));
Z{1}(i,j)=Theorem3(9.64,0.01643,0.163104,9,y1(i),9-y2(j));
Y{2}(i,j)=Theorem4_1(1817.2,0.53512,0.30758,1800,x1(i),x1(j));
Z{2}(i,j)=Theorem4_1(9.64,0.01643,0.163104,9,y1(i),y1(j));
Y{3}(i,j)=Theorem4_2(1817.2,0.53512,0.30758,1800,x1(i),1800-x2(j));
Z{3}(i,j)=Theorem4_2(9.64,0.01643,0.163104,9,y1(i),9-y2(j));
Y{4}(i,j)=Theorem4_3(1817.2,0.53512,0.30758,1800,x1(i),1800-x2(j));
Z{4}(i,j)=Theorem4_3(9.64,0.01643,0.163104,9,y1(i),9-y2(j));
    end
end
% meshc(x1,x2,Y{1})
% xlabel('sup_{[0,La)}U(t)'),ylabel('inf_{[0,La)}U(t)'),zlabel('Probability Density');
% title('The joint distribution of sup_{[0,La)}U(t) and inf_{[0,La)}U(t)');
% meshc(x1,x1,Y{2})
% xlabel('sup_{[0,Ta)}U(t)'),ylabel('sup_{[Ta,La)}U(t)'),zlabel('Probability Density');
% title('The joint distribution of sup_{[0,Ta)}U(t) and sup_{[Ta,La)}U(t)');
% meshc(x1,x2,Y{3})
% xlabel('sup_{[0,Ta)}U(t)'),ylabel('inf_{[0,La)}U(t)'),zlabel('Probability Density');
% title('The joint distribution of sup_{[0,Ta)}U(t) and inf_{[0,La)}U(t)');
% meshc(x1,x2,Y{4})
% xlabel('sup_{[Ta,La)}U(t)'),ylabel('inf_{[0,La)}U(t)'),zlabel('Probability Density');
% title('The joint distribution of sup_{[Ta,La)}U(t) and inf_{[0,La)}U(t)');

% meshc(y1,y2,Z{1})
% xlabel('sup_{[0,La)}U(t)'),ylabel('inf_{[0,La)}U(t)'),zlabel('Probability Density');
% title('The joint distribution of sup_{[0,La)}U(t) and inf_{[0,La)}U(t)');
% meshc(y1,y1,Z{2})
% xlabel('sup_{[0,Ta)}U(t)'),ylabel('sup_{[Ta,La)}U(t)'),zlabel('Probability Density');
% title('The joint distribution of sup_{[0,Ta)}U(t) and sup_{[Ta,La)}U(t)');
% meshc(y1,y2,Z{3})
% xlabel('sup_{[0,Ta)}U(t)'),ylabel('inf_{[0,La)}U(t)'),zlabel('Probability Density');
% title('The joint distribution of sup_{[0,Ta)}U(t) and inf_{[0,La)}U(t)');
meshc(y1,y2,Z{4})
xlabel('sup_{[Ta,La)}U(t)'),ylabel('inf_{[0,La)}U(t)'),zlabel('Probability Density');
title('The joint distribution of sup_{[Ta,La)}U(t) and inf_{[0,La)}U(t)');